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Abstract. We study the interplay of Anderson localization and interaction in a two chain Hubbard ladder 
allowing for arbitrary ratio of disorder strength to interchain coupling. We obtain three different types of 
spin gapped localized phases depending on the strength of disorder: a pinned 4&f Charge Density Wave 
(CDW) for weak disorder, a pinned 2kp CDW^ for intermediate disorder and two independently pinned 
single chain 2/cf CDW for strong disorder. Confinement of electrons can be obtained as a result of strong 
disorder or strong attraction. We give the full phase diagram as a function of disorder, interaction strength 
and interchain hopping. We also study the influence of interchain hopping on localization length and show 
that localization is enhanced by a small interchain hopping but suppressed by a large interchain hopping. 

PACS. 71.10Pm Fermions in reduced dimensions (anyons, composite fermions, Luttinger liquid, etc.) - 
71. 23. An Theories and models; localized states 



1 Introduction 

In the recent years, the two-chain Hubbard model received 
a lot of theoretical attention both analytically [ M, j[B 
10,110,0 and numerically Iff. 
The outcome of these studies was that the two-chain Hub- 
bard model with repulsive interactions would present a 
metallic phase with a spin gap containing dominant "d- 
wave" superconducting fluctuations, a situation very remi- 
niscent of the situation in two-dimensional cuprate super- 
conductors. This result leads to the conjecture that the 
■ two-chain Hubbard model could be considered as a toy- 
model of cuprate superconductivity. More recently, exper- 
iments revealed that the material with two-chain structure 
Sri4Cu2404i doped with Ca could become metallic pfjf 
and displayed a superconducting phase under pressure. 
This superconducting phase has been interpreted as the 
stabilization of "d-wave" fluctuations by interladder Joseph- 
son coupling. Such an interpretation is not clear enough 
because doping by Ca is not only introducing carriers in 
the system, but also inducing disorder. 

In a one-dimensional system, even an infinitesimal dis- 
order induces Anderson localization and an insulating state 
in the presence of interactions M]. Therefore it is crucial 
to understand the effect of interchain coupling on An- 
derson localization. At present, the effect of an infinitesi- 
mal disorder in the limit of strong interchain coupling has 
been analyzed |23,E3|. The outcome of this study is that 
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although the presence of the spin-gap suppresses the cou- 
pling of disorder to 2kp fluctuations, disorder can still cou- 
ple to Akp fluctuations (hp being the Fermi wave vector) 
leading to an insulating phase rather than a superconduct- 
ing phase at low temperature. Indeed, there is some ex- 
perimental evidence from transport measurements [^4] , |25| , 
^6|,^ that the phase at low temperature under ambient 
pressure is a pinned Charge Density Wave (CDW) state. 
This, as well as the disappearance of spin gap in NMR 
and the reduction of anisotropy of electrical conductivity 
close to the superconducting transition has lead to the idea 
that pressure could cause a dimensional crossover from 
ID insulator to 2D superconductor and that supercon- 
ductivity in this system is intrinsically two-dimensional 
phenomenon 1 28 , 2^,^, [u]. Building on the results of [p3[ , 
a phcnomcnological theory of transport in the insulating 
state of the ladder system was proposed in . 

Besides its relevance to ladder compounds, the study of 
two-chain systems is also important for the study of trans- 
port in systems such as metallic carbon nanotubes [ j33| , |34| ] 
and quantum wires |35, 36, 3l],|38|]. It has also a theoreti- 
cal interest in itself since the interplay of localization and 
interaction in dimension higher than one remains poorly 
understood) 39 and coupled quasi-one dimensional system 
may allow for a controlled study of the effect of increasing 
space dimensionality on disorder and interaction. 

In the present paper, we extend the work of Refs. [ p2| 
p3| to arbitrary ratio of disorder to interchain coupling. A 
study in the case of spinless fermions revealed an interest- 
ing transition from two chain to single chain behavior as 
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disorder strength was increased p0|| , indicating confine- 
ment of fermions by disorder. We shall analyze the two 
chain system with disorder by bosonization |ftl[|4^ .[l3|. [i~4| 
[45|| and RG (renormalization group) techniques and dis- 
cuss the different disordered phases that are obtained as 
well as the confinement phase |4(||^]]. We use the same 
techniques as Kawakami and Fujimoto (48), who analyzed 
the interplay of disorder and Umklapp scattering in the 
two chain ladder at commensurate filling in the regime in 
which interchain hopping is dominant p8| . However we ex- 
amine a different system at incommensurate filling, but al- 
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low disorder to become stronger than interchain hopping. 

Other related work in the ladder considered the effect of where C dw° 0) ( C dw« (x)) denotes a operator of CDW 



forward scattering on the renormalization of disorder |}49| , 
neglecting the influence of gap openings on local- 
ization properties. Such an approximation is restricted to 
the limit of high temperature or short system, and can be 
recovered from the general RGE (renormalization group 
equation) of the present article. 

The plan of the paper is as follows. In Sec. || we will 
recall the derivation of the bosonized Hamiltonian of the 



(charge density wave) with in phase (out of phase) be- 
tween two chains and is defined by 



Ocdw°( x ) = e l * p+ {sin0 p _ cos <t> a+ cos cp 
— i cos 4> p - sin (f> a+ sin } , 
Ocdwt (x) — e^ p+ {cos 9p- cos a - cos (j>, 
+ i sin Op- sin a - sin <j) a+ } . 



(7) 
(8) 

system. Then, in Sec. |we will show the RG equations T , , tt t \ ■ ■ ui n 

J . ,, i -i ,i c - 1 r ,i i- •; r i ■ , i • -Lhe quantity ll ur (x) is a variable canomcally conjugate 

of the system while those of the limit of large interchain 

hopping is derived in Sec. 4 to compare with other re- 
lated work. In Sec. [|, we will demonstrate several types of 
spin gapped phases based on the outcome of RG. Finally 
discussion is given in Sec. ^. 
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2 Hamiltonian 

The original lattice Hamiltonian of two-chain Hubbard 
model with disorder reads: 
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From bosonization rules, these quantities are written as 

(to = 0, 71") 
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We turn to the bonding antibonding band representa- 
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and use the continuum limit to apply bosonization techniques [| 
23Lp0fl. Thus the Hamiltonian can be rewritten as: 
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The coupling constants of the bosonized Hamiltonian 
(^) can be expressed as the function of those of the lattice 
Hamiltonian (pi). They are given by 
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K v - = 1, 

U v - = Up, 
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and 

(21) 

(22) 
(23) 
(24) 
(25) 
(26) 
2Dld(x - 

(x') = 2D b r S(x - x')S r , r >- v F is the 
Fermi velocity in the absence of interaction. We note that 
there is a freedom of sign in choosing the initial condition [f34 
although the result remains the same. Actually, the initial 
value of yu y§ and yio in Eq. ( p(j| ) has a sign opposite to 
that of p7| , i.e, the correspondence is given by the re- 
placement of </> CT + — > (j> a + + ir/2. We consider only the 
backward scattering for impurity, which is relevant to the 
localization. 



3 RG equations 

The RG equations can be derived from the Hamiltonian 

using the Operator Product Expansion (OPE) approach] 
In the RGE, we neglect the velocity differences between 
the different sectors. To be consistent with such an ap- 
proximation, we also have to neglect 2/4, ys in the RGE 
Failing to do so leads to unphysical results such as a spin 
gap for ij_ — and no disorder in disagreement with the 
exact results on the single chain Hubbard model|Q. It 
is also necessary to expand the RG equations of K a ± to 
first order to ensure SU(2) symmetry. The renormalization 
group equations read: 
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where we have introduced y a ± and i± defined respectively 
by K a± = 1 + y a± /2 and = ^7. 

In these RG equations, we have neglected the genera- 
tion of a 4k f disorder from 2k f disorder since the gener- 
ation of such term in RG demands special treatment p3[ . 
In the present calculation, we will be finding SC d and SC S 
phases, in which the 2kp disorder is regarded as irrelevant. 
However, based on the study of Ref. [E3[ we already know 
that 4k p disorder is always relevant in the SC d phase and 
relevant in the SC S except for K p+ > 3/2. Therefore, in 
reality these phases of SC d and SC S will be replaced by a 
pinned 4kp CDW as discussed in Ref. p8~|. 



3.1 The limit of V = 0: the pure two-chain ladder 

In the limit where V = 0, a simplified form of the RGEs 
is obtained as: 
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Kp- {(vl +vl + y7)Jo(4:K p -t±)}] , (43) for small I, one should have y a - < 0. This in turn makes 
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These equations are identical to those previously de- 
rived in . Taking the limit of t± — > +00, which is the 
limit of the two chain system, these equations are further 
simplified(l|,0: 
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Obviously, one has: A p _ — > +oo,j/ CT+ — » —00. More- 
over, it is easily seen from initial conditions that originally 



1 2/2 1 < 1 2/3 1 and |y 9 | = \y w \ , resulting in 
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2/2,2/10 less relevant and results in both — -w 
coming more negative. As a result, one hnds — ► —00. 
Thus p _ and become locked in the ladder system. 
Further analysis of these RGE shows that there are two 
possible fixed points [^4[01 depending on the sign of U 
as represented on table [l]. The fixed point for U > 
correspond to the formation of d-wave superconductivity, 
whereas the one for U < corresponds to the formation 
of s-wave superconductivity. 



3.2 The case t± = 0: single chain limit 

This case corresponds to two decoupled disordered Hub- 
bard chains and has been analyzed in details in Ref . [^IJ . A 
crucial test of the validity of our RG equations is whether 
we can recover the RG equations of Ref. plj by putting 
t± = in Eqs. For t± = 0, it is straightforward 

to show that the solution of the RG equations (|27j)-(41 ) 
satisfies: 
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As a result, the RG equations (^)-(^T)) are reduced to the 
following equivalent set: 
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Using the expansion: A p+ = 1 + y p +/2, we rewrite the 
RG equation for A p+ in the form: 



dy P 



dl 
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with y p +(0) = —?/(0). It is then easily seen that yi(l) = 
(?/p+(0 +2/(0)/ 2 , f or all I. This allows to rewrite the RGE 
for T> as: 



1 dV „ r ^ 
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The RGE equations are thus reduced to those of the single 
chain case (see Eq. (3.12) in Ref. §l|). Then, our RG 
equations have the correct limit of t± = 0. Obviously, for 
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T> = 0, we recover the RGEs that can be applied to a single 
chain system. In particular, we can easily check that the 
system develops a spin gap for U < but remains gapless 
for U > 0. Thus, we see that the RGEs ©-© correctly 
capture the limit of a single chain system. 



4 The limit t± — > oo: two chain regime 

In this section, we consider the limit t±_ — > oo of the equa- 
tions (|27 )-(|4ll) in order to make contact with the work of 
Refs . fZlfc 5 1|, [49|, jsp) . In the limit t± — > 00, the Bessel func- 
tions are negligible, so that Eq. (p7|)-(pl|) are reduced to: 



disorder and large t± thus read: 
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These equations imply that antisymmetric disorder is more 
relevant than symmetric disorder in the case of attrac- 
tive interaction and less relevant for repulsive interactions. 
Such result can be understood in the following way: for re- 
pulsive interactions, the system tends to form interchain 
Cooper pairs and thus to have a repartition of charges 
that is symmetric between the two ladders. This makes 
the system less sensitive to antisymmetric disorder. For 
attractive interactions, the system prefers to form intra- 
chain Cooper pairs. If electrons of the Cooper pair can tun- 
nel on the opposite chain in a virtual process, the energy 
of the Cooper pair is lowered. This results in interchain 
repulsion of the pairs and thus antisymmetric CDW fluc- 
tuations. This mechanism obviously makes antisymmetric 
disorder a more relevant perturbation. 

Comparing with Ref. |48[ ], we have similar equations 
once we put the Umklapp process of Ref. 0] equal to zero. 
However, an important difference between the present pa- 
per and [Esj is that in |Q a random t± term is also in- 
cluded in the Hamiltonian. This random ij_ can result in 
closure of the transverse charge gap and thus can lead 
to a rather different physics from the one discussed here, 
closer to the one considered in We will com- 

pare in more details our results and those of [Q in the 
next section. 



5 Results 



From RG equations d^)-©, we can distinguish four dif- 
ferent phases where the corresponding fixed points are 
shown on table 2 ( explained below). 



dV s 
dl 
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(85) 

Let us compare our equations with those in 1 50 . Equa- 
tions (11), (12), (14), (15) in Ref. J5Q] correspond respec- 
tively to equations (gj), (|8^), (p3[ ) and j79| ) in the present 
paper. However, in Eqs. (12), (14) and (15) of Q we do 
not find terms arising from coupling between interactions, 
which give rise to the same order. More importantly, Eq. 
( p5| ) is in disagreement with Eq. (10) in Ref. [|50| . The rea- 
son is that in Ref. | p0fl , the 2/5 term was kept although it 
should drop from the Hamiltonian when t± becomes suf- 
ficiently large. The correct initial RG equations for weak 



5.1 SCd Phase 

For large repulsive interaction and small disorder, we ob- 
tain the SCd phase, where strong coupling in interaction 
is reached before strong coupling in disorder. The result- 
ing flow of coupling constants is represented on Fig. 0. The 
ground state has (9 P -) = 0, {(/>„-) = (4>a+) = § (I in table 
2). This phase to the approximation of RG is a metallic 
phase. It has been discussed in Q as case A. However, 
including coupling to Akp disorder, it is possible to show 
that such a phase is in fact localized |2^,^3|, 48 1 due to 
pinning of 4fcp CDW fluctuations by disorder [22,p3|. We 
can see that in this phase, D s is more relevant than D a 
which is consistent with Eqs. (|86|)-(|87|). For stronger in- 
teraction, this behavior can also be understood from Eqs. 

■ Replacing p _, (j> a ± by their expectation values in 
these equations, we obtain Ocdw" ~ e l ^ p + cos <j) P - and 
Ocdw* ~ e l ^"+ (cos a - cos0 CT+ + isin6> p _ sin0 CT _) where 
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we have kept only the operators with power-law or expo- 
nential decay. Clearly, Ocdw™ contains more disordered 
operators with exponential decay than Ocdw/o. As a re- 
sult, CDW° are dominating over the CDW^ in the SCd 
phase, in agreement with the picture of interchain pairing. 
This also implies that the random potential that couples 
to Ocdw is more relevant than the one that coupled to 
Ocdw" even deep in the SCd phase. 



5.2 SCs phase 

For interaction being moderately attractive and disorder 
being small, we obtain the phase where strong coupling 
in interaction is reached before strong coupling in dis- 
order. The resulting flow of coupling constants is rep- 
resented on Fig. ^. In this case, the ground state has 
(0 P _) = {<j> a -) = (<j> a+ ) = (lib in table 2). Similarly 
to the SCd phase, the SCs phase is localized by disorder 
due to the presence of subdominant 4/cp fluctuations [22, 
|2l| . However, for attractive enough interactions ]22|, 23 1 
the Akp CDW can be depinned by quantum fluctuations 
in contrast with the SCd phase. It can be seen that in the 
SCs phase, D a is increasing faster than D s in agreement 
with Eqs. (|86|)~ (|87|). The persistence of this behavior even 
when the spin gap is well f ormed can be understood by the 
same arguments as in 5.1. Comparing Eqs. (Q) and (ph. In 
Eq. (||), only the operator cos6> CT _ has exponentially de- 
caying correlation, whereas in Eq. (j7j) the three operators 
cosd 



sin</> CT _ and sin</> ff + are disordered with exponen- 
tially decaying correlations. As a result, the CDW have 
a much faster exponential decay than the CDW 7 ' fluctu- 
ations inside the SCs phase. In turn, this makes D s much 
less relevant than D a . Physically, this can be understood 
by noting that attractive interactions form Cooper pair 
inside a single chain, leading to subdominant CDW 7 ' fluc- 
tuations whose transverse charge fluctuation is compatible 
with that of SCs fluctuations. 



5.3 PCDW^ phase 



5.4 Confinement 

If strong coupling in disorder is reached before t±(l) = 
1, there is a possibility of confinement of carriers in the 
transverse direction corresponding to the irrelevance of 
t± (IV in table 2). This is represented on Fig. ^ 

(curve (7))- Note that on the figures t± stands for t±. 
This confinement can be understood as the formation of 
a localized state at an energy scale higher than t±. This 
time, a spin gap is formed at an energy scale higher than 
ii. As a result, we have to distinguish two types of pinned 
2k F CDWs. One is the pinned 2k F CDW of the two 
chain system (III in table 2), and the second one is the 
combination of two pinned 2k F CDW in each chain with 
irrelevant interchain hopping (IV in table 2). In the first 
case, there is a definite phase relation between the CDW 
in chain 1 and the CDW in chain 2, while in the second 
case, such relation is less definite. Such regime could not 
be obtained with the approximations used in [H or [^3| 
which are valid only when t± is relevant. Further, we note 
that confinement can also result from the presence of an 
attractive interaction in the absence of disorder as seen 
curves (1), (2) and (3) in Fig. |) (Ha in table 2). 

Here the difference in the disordered state between Ha 
and IV is examined besides the mechanism for confine- 
ment. In the region Ha of Fig. [|, a Josephson coupling 
is present and dominates the low energy properties. The 
physics of this region can be understood by describing the 
Cooper pair as hard core bosons. The problem is reduced 
to a disordered bosonic ladder with interchain hopping 
considered in Ref. |35). Thus it is found that the phase in 
region Ha is a pinned Ak F CDW. In region IV, the local- 
ization effects are stronger than the Josephson coupling, 
leading to two independently pinned 2k F CDWs. 



5.5 Criterion for transitions between the different 
phases 

We examine boundaries between several types of phases 
, which are shown in Fig. [5[ 



For intermediate disorder, we obtain the PCDW" phase 
where strong coupling in the RG equations is obtained 
for disorder first. The flow of the coupling constants is 
represented on Fig. [| In this phase, one finds that D a 
is diverging faster than D s . This corresponds to the long 
range ordering with (8 p -) = {<j) a +) = {0 a -) = (III in 
table 2). This phase is the pinned CDW 7 ' of Ref. [E3[ and 
has also been discussed within a RG approach in p8fl . We 
note that such a phase appears only for disorder strong 
enough to compete with interaction ( i.e., not in the limit 
of infinitesimal disorder) for the present model. The pres- 
ence of such phase in the ladder system can have impor- 
tant consequences for transport properties since pinning 
properties of the CDW 7 ' are very different from the Ak F 
CDW. 



5.5.1 Confinement transition 

A criterion for disorder confinement can be obtained by 
comparing the length scale k oc . ( V(li oc ) — 1) with the 
length scale l 2c h. ( t±(l 2c h.) = !)• Clearly, if k oc , < l 2c h. 
the system is localized before entering into the two-chain 
regime and thus the ground state will be formed of two 
chains pinned on their own disorder. For weak interac- 
tions, this criterion corresponds to T)(0) > ij_(0). As can 
be seen on fig. ^, it gives the correct boundary between the 
pinned 2k F CDW (phase (III)) and the independently 
pinned 2k F CDWs (phase (IV)). For stronger interactions, 
the criterion for confinement is Ei oc . = (v F /a)e~ l,oc - = t±. 
Since the pinning energy is enhanced in the presence of 
interactions, this implies that interactions enhance con- 
finements effects as can be seen on Fig. |[ In the case of 
confinement by intrachain interactions one has to compare 
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the scale at which y(l) — — oo (for T> = 0) to the scale at 
which t±(l) ~ 1. Since we have: 



2/(0 



2/(0) 



l + y(0)Z' 



we find that confinement occurs when: 



2/(0) = 



ln(t x (0)) 



(88) 



(89) 



For i± = 10~ 2 , this gives confinement when Ua/nvp = 
—0.22, which is the correct value for Fig. when V < 

io- 8 . 

In Fig. [|, we see that the two confinement regimes 
merge to give rise to a single confined phase. An important 
question is whether the transition between the confined 
and the deconfined phase is a true phase transition (i. c. 
whether one can find an order parameter for the confined 
or deconfined phase) or if it is only a crossover. In the 
case of a commensurate potential ( at half-filling), it is 
well known |47],|56|,^7j that no order parameter exists than 
can distinguish the confined from the deconfined phase. 
Thus, wc do not expect to see a true transition in the 
disordered phase. Instead, we expect a progressive loss of 
phase coherence between the two CDWs in chain 1 and 2. 



region (I) and (III) is also obtained from the condition 
?/i (oo) = since the fixed point for y\ is given by oo for 
region (I) and — oo for region (III). The boundary between 
(III) and (lib) ( or between (IV) and (lib)) is also given 
by 2/2(°o)/y 3 (oo) = 1 since 2/2(00) / '2/3(00) > 1(< 1) for re- 
gion(III) and (IV) (region (lib) and (Ha)). The boundary 
between the confined and deconfined phase is given by the 
condition 



tl(/o) = I- 



(91) 



The quantity l c denotes a value at which ti takes a maxi- 
mum. The boundary between confinement and deconfine- 
ment in the pinned CDW region moves continuously to 
that in th SCs state. The former is determined by disorder 
while the latter is determined by attractive interaction. 
Such a continuous variation of the boundary indicates a 
fact that the confinement occurs by the combined effect of 
interchain hopping, intrachain interaction and disorder. 

Another global phase diagram is shown in Fig. ^ in the 
(t_L, U) plane with a fixed D = 1CP 5 where the notations 
are the same as Fig. 0. The effect of interchain hopping 
on confinement is opposite to that of disorder in Fig. 
The t± dependence on the boundary between (III) and 
(Ha) ((III) and (lib)) is small while it is noticeable for the 
boundary between (III) and (IV) ( (III) and (I)). 



5.5.2 Transition between PCDW" and SC phases 

This transition is obtained in the deconfined regime. A 
criterion for this transition can be obtained by comparing 
the 2k p CDW pinning energy with the gap in a— sec- 
tor. When the pinning energy is above the gap, the sys- 
tem gains more energy by being in the PCDW phase. The 
transverse gap in a— is obtained by applying single chain 
RG for t±e l <C 1 and two chain RG for t±e L ^S> 1. Com- 
paring the transverse gap with the pinning energy leads 
to a good agreement with the phase boundary between (I) 
and (III). 

5.6 Global phase diagram 

The global phase diagram is shown in Fig. || with the fixed 
t± = 10~ 2 . Based on the criterion in the previous subsec- 
tion, the boundary between pinned CDW phase (region 
(III)) and superconducting phases (region (I) or region 
(Ha)) is given by 

tf«r-(Jloc.) = 0. (90) 

The quantity li oc . is the scale at which one of the cou- 
pling constants of disorder is equal to one. This is also 
understood from Figs. |l|, ||, and || where the fixed point 
of y a - is given by —00 for regions (I) and (Ha) and 00 
for regions (III). In the confined region, the corresponding 
boundary is given by y a -(oo) = 0. With increasing D, the 
pinned CDW region is enlarged. The region for SCs state 
is much larger than that for SCd state since the spin gap 
exists even for the single chain. The boundary between 



6 Discussion 

We have examined the effect of interchain hopping on dis- 
ordered two-coupled chains using RG equations that en- 
able us to describe the crossover from the single chain to 
two chain regime. Using the RG equations, we have shown 
that a 2kp CDW 71 " could be obtained when disorder was 
strong enough compared to interaction. 

In the case of two chains of spinless fermionsjfO), it is 
known that for repulsive interactions there is a strong rein- 
forcement of Anderson localization, but complete delocal- 
ization for attractive interactions. By contrast, in the case 
of spinful fermions, the localization-delocalization transi- 
tion is known to be obtained only for K p+ > 3/2 i. e. large 
attractive interaction strength|23J. Moreover, in the case 
of repulsive interactions, in contrast with the spinless case, 
there is no enhancement of localization for infinitesimal 
disorder[^3| in the spinful system. This behavior is due to 
the fact that both for repulsive and attractive interactions, 
the localized phase in the presence of infinitesimal disor- 
der is a pinned 4fc^ CDW. Our RG study shows that when 
disorder becomes stronger (see figure ||) or interchain hop- 
ping becomes weaker (see figure ft), a 2kp pinned CDW" 
can be induced in the place of the AkF CDW. This pinned 
CDW" has a much shorter localization length than the 
pinned AkF CDW. Interestingly, the domain of existence 
of the PCDW" is larger for repulsive interactions. Thus, 
for weak but not infinitesimal disorder, localization will 
be reinforced by repulsive interactions. A similar effect is 
also obtained by reducing t± . This results from the renor- 
malization of interactions in the a— sector by disorder, 
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and could not be obtained within the infinitesimal disor- 
der limit of [p3f . Such behavior is represented on Fig. [7] 
where the localization energy i?i oc . = ■^ E - and the gap 



energy E gap 



_ V F 
£gap 



(6< 



and £ gap = ae gap being 



^loc. ^t- oxikj. s,gap 

respectively the lengthscale at which disorder or interac- 
tion become of order one) are represented a a function of 
tj_ for fixed disorder and interaction. A slight enhance- 
ment of the pinning energy can be observed for t± large 
enough to deconfine but not strong enough to make the 
gaps robust to disorder. 

Let us now turn to the question of the crossover from 
single chain to two chain behavior. In the limit of t± = 
i. e. the single chain fermion system, it has been known 
for some timep8| that the localization length was £i oc . ~ 
(1/ 'D) 1 /( 3 ~- R "p) for attractive interactions, but 
60c. ~ (l/D) 1 ^ 2 "-^' for repulsive interactions, leading 
to a reinforcement of localization by attractive interac- 
tions. This reinforcement of localization results in a max- 
imum of localization length in the vicinity of the non- 
interacting point. For large t± and infinitesimal disorder, 
it has been shown in p3[ that the localization length was 
£/oc. ~ (l/D) 2 /( 3 ~ 2Kp +> so that the reinforcement of lo- 
calization by attractive interactions disappeared. In the 
present case of intermediate t± and non-infinitesimal dis- 
order, the behavior of the localization length deduced from 
the RG is richer. For small interaction strength, disor- 
der is strong enough to prevent the formation of SCs or 
SCd phase in the system, and the behavior of localization 
length is essentially one dimensional, leading to a max- 
imum of localization length in the vicinity of the non- 
interacting point. For stronger interactions, the SCs and 
SCd phases can develop, leading to a strong enhancement 
of localization length caused by the weaker pinning of Akp 
CDWs compared to 2kp CDWs. This result in the pres- 
ence of two minimums of the localization length as a func- 
tion of interaction. The corresponding behavior of the pin- 
ning energy E\ oc . — «f/£ioc. in the case of D = 10 -5 and 
Ua/(nvp) — 0.5 is represented on Fig. || .Interestingly, 
the presence of these minimums is asymmetric especially 
for small t±. This behavior can be understood by looking 
in more details at the RG equations. At a scale I < ?2cft.j 
where l 2c h. is such that t±(l2 C h.) = ^f-, the RG equations 
are those of the one dimensional system, and interactions 
in the spin sector are renormalized towards zero in the 
case of repulsive interactions. Only for I > l^ch. are inter- 
actions in the spin sector growing again to lead to the spin 
gap. By contrast, for attractive interactions, the interac- 
tions are always growing in the spin sector, irrespective of 
the strength of t ± . This implies that the spin gap is much 
weaker in the case of small ij_ and repulsive interactions 
thus making the system more sensitive to disorder. There- 
fore, in the case of the ladder, enhancement of localization 
is not determined only by the sign of interaction but also 
by their strength. In the case of weak interactions, the 
behavior is similar to the single chain case. For stronger 
interactions, the effect can be inverted and attractive in- 
teractions can have a delocalizing effect. Finally, for even 
stronger interaction, the difference between attractive and 
repulsive interactions becomes unobservable. It would be 



interesting to study how the charge stiffness is modified 
by disorder in these different regimes. 

A problem related to the single chain to two chain 
crossover is the issue of confinement i.e. the irrelevance of 
single particle hopping under RG. As discussed in section 



[5.5.1 , the criterion for deconfincmcnt is -E10C. = t±. This 
criterion indeed reproduces the phase boundary between 
the phase III (deconfined) and IV (confined) as can be 
seen by comparing Figs. || and pi Interestingly, a kind 
of confinement seems to be observed in phase (Ha) as a 
result of attractive interaction. However, an effective two- 
particle hopping is known to be generated in this case |39[ 
p0{ and plays the role of a Josephson coupling between 
the chain. Such coupling causes the formation of a gap in 
the transverse charge mode so that the phase (Ha) is still 
the SCs phase in the absence of disorder and becomes the 
4fcp pinned CDW in the presence of disorder. Confinement 
of carriers by disorder is obtained in the case of strong 
disorder. Such a confinement has not been considered in 
the previous study of the spinless fermions casejio] due to 
a choice of large t±. The difference in confinement is as 
follows. In the spinful case and for a single chain, disorder 
generates a spin gap in the localized phase. Thus, in the 
localized phase (region (IV) in Fig. 5), a gap in the single 
chain has to be overcome to permit single particle hopping. 
If one turns on a weak single particle hopping between two 
chains, its effect is negligible since the fields <p p+ in both 
chains are disordered as a result of the presence of strong 
disorder (this is to be contrasted with the phase Ha). In 
the case of spinless fermions, no spin degree of freedom is 
available to induce a spin gap so that the "confinement" 
results only from the growth of the yj_ term under RG. As 
one can see from this discussion, this effect is not limited to 
the two chain problem but is also present for any number 
of coupled chains. Therefore, confinement by disorder is 
a generic feature of coupled disordered spinful chains and 
should also be observable for instance with three coupled 
chains. 



7 Conclusion 

In the present paper, we have applied RG techniques to 
study Anderson localization in a two-chain Hubbard lad- 
der for arbitrary ratio of t±,U,T>. We have found that 
there were three different types of localized phases: the 
Ak F pinned CDW (I and II) the 2k F CDW 77 (III) and the 
decoupled single chain 2kp CDW (IV). We have been able 
to obtain the full phase diagram of this system. We have 
shown that a weak interchain hopping could lead to a rein- 
forcement of localization, but a strong interchain hopping 
resulted in a reduction of localization as discussed in | p3| . 
The presence of 2kF CDW phases should have interesting 
consequences on transport properties that need separate 
investigation. A remaining open problem is to study the 
generation of the Akp disorder (i.e. disorder that couples 
to the Akp CDW operator) directly from the RG equa- 
tions and the resulting overall behavior of conductivity. 
This will be left for a future study. 
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Table 1. The two possible fixed points in a two-chain Hubbard 
ladder 



U > 


<*,-> = o 


<^-> = f 


<^+) = f 


U < 


(Op-) = o 







Table 2. The fixed points of the disordered two-chain Hubbard 
ladder 



I 


(Op-) = o 


(<t>*-) = f 




fx —* 00 


Ha 


K p - ~ 1 


K(j— ^ 1 


(<^+> = o 


fx -> 


lib 


<^-> = o 


=0 


(<?W) = o 


fx — > oo 


III 


(Op-) = o 




<<^+> = o 


fx — » oo 


IV 


AT„- =± 1 


K a - ~ 1 


(4>*+) = o 


fx -> 



-1 



.iii! 












■ y c - 

■ . , , , i 


A, 



0.02 
D 
0.01 



-0 



/ 



10 



Fig. 1. RG flow as a function of I for the SCd phase with 
Ua/irv-p = 0.5, fx — 10~ 2 and D — 10~ 7 where the left axis is 
for j/i and 3/ CT _ and the right axis is for D a (dashed curve) and 
D B (dotted curve). 




Fig. 2. RG flow as a function of I for the SCs phase with 
Ua/nvF = —0.1, fx = 10~ 2 and D — 10~ 7 where the notations 
are the same as Fig. 1. 




D 



Fig. 3. RG flow as a function of { for the the PCDW phase 



with Ua/-KV F = 0.01, fx = 10~- 
notations are the same as Fig. 1. 



and D = 10 _c> where the 




/ 



Fig. 4. RG flow of fx showing both confinement (irrelevant 
fx) and deconfinement (relevant fx) with Ua/nvF = -0.5(1), 
-0.4(2), -0.3(3), -0.2(4), -0.1(5) and 0(6) where t± = 10~ 2 
and D — 10~ 3 . The curve (7) also shows the confinement for 
Ua/iTv F = 0.1, fx = 10~ 2 and D = 10" 1 . 
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Fig. 5. Phase diagram in the (log(X>),{7) plane for a fixed 
t± — 10~ 2 . The respective phases denote the SCd state ( region 
(I)), the confined SCs state (region (Ila)), the SCs state (region 
(lib)) the pinned CDW state (region (III)) and the pinned 
CDW state with confinement (region (IV)). 
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Fig. 7. Behavior of 2Uf pinning energy (inverse of localization 
length) and gap energy as a function of interchain hopping for 
y = 0.5 and V = 1CP 6 . For a given value of t± only the largest 
of the two energies is represented. We note that for small in- 
terchain hopping, there is a small reinforcement of localization 
with respect to the single chain limit. The origin of this re- 
inforcement is the development of 2kF CDW^ fluctuations in 
the system. With stronger interchain hopping, derealization 
effects become apparent. Finally, when the interchain spin gap 
becomes larger than the 2&f pinning energy, the system devel- 
ops a 4k f pinned CDW with a much smaller pinning energy. 
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Fig. 6. Phase diagram in the (t±,U) plane for a fixed D = 
1CP 5 where the notations are the same as Fig. 5. 
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Fig. 8. The behavior of the pinning energy (inverse localization 
length) a a function of interaction strength for D = 10 -5 and 
t± = (single chain case, solid line) and t± =0.1 (two chain 
case, dotted line). For weak interaction, the pinning energy is 
the same for both systems. This corresponds to the confine- 
ment regime. For stronger interaction, the two chain system is 
always less localized than its single chain counterpart. Thus, 
a large enough t± always induces derealization. For an even 
stronger interaction, the two chain system goes into the pinned 
AUf CDW phase with a much lon ger pinning length. This limit 
has been previously discussed in p3[. 



